State dependent effective interaction for the hyperspherical 
formalism with noncentral forces 



Nir Barnea 1 , Winfried Leidemann 2 ' 3 , and Giuseppina Orlandini 2,3 
1 The Racah Institute of Physics, The Hebrew University, 91904, Jerusalem, Israel 
2 Dipartimento di Fisica, Universitd di Trento, 1-38050 Povo, Italy; 
s Istituto Nazionale di Fisica Nucleare, Gruppo collegato di Trento. 

(February 8, 2008) 



Abstract 

The recently developed effective interaction method for the hyperspherical 
harmonic formalism is extended to noncentral forces. Binding energies and 
radii of three- and four-body nuclei are calculated with AV6 and AV14 NN 
potentials. Excellent results for the convergence of the expansion are found, 
particularly for the three-nucleon system. Due to the higher density the con- 
vergence rate is a bit slower for the alpha particle. In comparison to central 
potential models there is only a very slight deterioration of the convergence 
due the tensor force, while other potential terms have no visible effect on the 
convergence. The obtained values for binding energy and radii also agree well 
with the results in the literature obtained with other few-body techniques. 
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I. INTRODUCTION 



In the last years great progresses have been achieved in microscopic calculations of few- 
body systems. Quite a few different exact approaches have been developed. Among those the 
effective interaction method is particularly interesting since it is a typical many-body rather 
than a few-body approach. In fact effective interactions are used in shell-model calculations 
of complex nuclei, where the wave functions are expanded in a model space consisting in a 
truncated single particle harmonic oscillator (HO) basis. The introduction of an effective 
rather than a bare interaction is necessary in order to compensate for the truncation of the 
space. A few years ago such an effective interaction method has also been introduced in 
few-body physics. So called no-core shell-model calculations have been carried out, where 
one keeps all the nucleons active and where one also removes the spurious center of mass 
motion finally leading to HO basis functions that depend on the Jacobi coordinates I0J. 

Recently we have proposed a different effective interaction method |J using instead 
of the HO basis an expansion of the wave function in hyperspherical harmonics (HH). In 
the HH formalism the Jacobi coordinates are replaced by a single length coordinate, the 
hyperradius, and a set of 3A — 4 hyperangles. The HH are the A-body generalization of 
the 2-body spherical harmonics, and likewise depend only on the hyperangular (angular) 
coordinates in the hyperspherical (spherical) decomposition of the A-body (2-body) system. 
In general, the wave function can be expanded in a series consisting of products of HH basis 
functions and hyperradial basis functions. The HH basis is widely used in the calculation 
of few-body wave functions, though the convergence is very often problematic. In order to 
improve the convergence one may introduce proper correlation functions (see, e.g., [|H^]). 
The correlation functions, however, lead to various undesirable features in the calculation. 
One not only looses the orthogonality of the basis functions but also has no unique way 
to determine the correlation functions. Therefore, as an alternative approach, we have 
studied the HH effective interaction (EIHH) method to accelerate the slow convergence of 
the expansion. Employing simple central NN potential models as bare NN interaction for 
light nuclei in the mass range of A = 3 — 6 we could show that the EIHH method leads to 
an extremely rapid convergence for the investigated observables. One important reason for 
the excellent convergence behavior is due to the fact that the effective interaction becomes 
state dependent in a natural way, i.e. the interaction does not only depend on the state of 
the two interacting particles, but also on the state of the rest system. 

Due to the excellent convergence the EIHH method has proven to be a very promising 
tool in few-body physics. On the other hand, in Ref. |2| only central potential models have 
been considered while realistic NN interactions are much more complicated and in principle 
could deteriorate the convergence rate. Therefore it is of great importance to check whether 
the EIHH method leads to similarly good results also for more realistic NN potentials. In 
this paper we investigate the method for two noncentral NN interaction models (AV6, AV14 
0) calculating ground state energies and radii for the A = 3,4 nuclei. To this end the 
convergence patterns are studied in detail and the final results are compared to results in 
the literature. 

The paper is organized as follows. The EIHH method for noncentral forces is described in 
Sec. [TI|. Numerical results, including a comparison with results obtained with other few-body 
techniques, are given in Sec. ITTTI and conclusions are drawn in Sec. [TV. 
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II. THE HH EFFECTIVE INTERACTION 




Before introducing the effective interaction we first give a short summary of the HH 
formalism. To illustrate the hyperspherical coordinates we start from the definition of the 
center-of-mass coordinate R — \ J2f=i U an d the normalized reversed order N = A — 1 
Jacobi coordinates 

/ A 'a i A 

(1) 

The Jacobi coordinate ffj consists of a radial coordinate rjj and a pair of angular coordinates 
rjj = (8j, (f)j). These coordinates are then transformed into the hyperangular coordinates 
a 2 , • • • j through the relations 

n 

sin a n = r) n / p n , p 2 n = p 2 n _ x + rf n = Vj ■ ( 2 ) 

3=1 

Note that the hyperradial coordinate p = p^ is symmetric with respect to permutations of 
the underlying single particle coordinates. The 3N = 3(A — 1) internal coordinates for the 
A-particle system consist of the hyperradial coordinate p and the 3iV — 1 "hyperangular" 
coordinates £l N = {fjx, fj 2 , • • • , fj N} a 2 , a 3 , • • • , a N }. 

With the hyperspherical coordinates one can write the Laplace operator for n Jacobi 
coordinates n = 1 . . . N, as a sum of two terms 

n ~ P^Opf" d Pn 7n n ' ^ 

The hyperspherical or grand angular momentum operator K\ of the n Jacobi coordinates 
can be expressed in terms of the squared angular momentum associated with the nth Jacobi 
coordinate, and K^-i as follows [§] 

d 2 3n - 6 - (3n - 2) cos(2a n ) d 1 >2 1 p. ... 

aa A n sm{2a n ) oa n cos^ a n sin a n 

where we define K\ = . The total angular momentum operator associated with n coordi- 
nates is L n = L n _i +l n - The operators K 2 , K 2 _ ly L 2 n and L nz commute with each other. 
The hyperspherical harmonic functions 3^[jr n i are the eigenfunctions of K 2 with eigenvalues 
K n (K n + 3n — 2). The symbol [if n ] stands for a set of quantum numbers including, e.g., 
£±, ...,£ n , L 2 , ■■■,L n and K 2 , ...,K n . For more detailed information on the HH functions see 
II- 

The HH functions do not possess any special properties under particle permutation. 
Therefore the first step in applying the expansion to the A-body problem is the symmetriza- 
tion of the HH basis. In the current work we employ two powerful algorithms |iO| , |TT| recently 



developed for the construction of a HH basis with well defined permutational symmetry. 

Employing noncentral potentials it is usually more convenient to work in a jj coupling 
scheme. However, constructing symmetrized basis functions in this coupling scheme is more 
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complicated, roughly by a factor 2 A , in comparison to the LS scheme. Therefore we work in 
the LS scheme coupling HH basis functions with definite permutational symmetry with spin- 
isospin basis functions into anti-symmetric many-body states with total angular momentum 
J. 

Now we turn to the question of the effective interaction. In general we would like to use 
the HH basis functions to solve the A-body Hamiltonian 

A -2 A 

# = £|b + E^, (5) 
i=i z '' n i<j 

where Pi is the momentum of the iih particle and m the nucleon mass, while 

V*i = V(uj) = £ V^ir^d^ihj ; ffi, 3,, n, r 3 ) (6) 
p 

is the NN interaction which also depends on the spin (<jj) and isospin (ji) operators of the 
two interacting particles. As already mentioned one usually introduces correlation functions 
in order to obtain sufficiently converging results In our alternative EIHH method of 

Ref. we use instead of the bare NN interaction an effective interaction inside the model 
space. Defining P as the projection operator onto the model space and Q = 1 — P as the 
projection onto the complementary space, the model space Hamiltonian can be written as 



Hr 



&2m 



P + P 



Kj 



P. (7) 

eff 



In the HH formalism the model space can be defined as a product of the hyperradial subspace 
and the set of HH basis functions with generalized angular momentum quantum number 
K < K 

In general the effective interaction appearing in Eq. (|7|) is an A-body interaction. If 
it is determined without any approximation, then the model-space Hamiltonian provides a 
set of eigenvalues which coincide with a subset of the eigenvalues of the original full-space 
Hamiltonian, Eq. (0). Since the effective interaction has to contain the whole A-body 
information of Q-space the calculation of this A-body effective interaction is as difficult 
as finding the full-space solution. Because of these difficulties we introduced in Ref. 
instead of the exact A-body an approximate two-body effective interaction which, however, 
satisfies the following properties: (i) Viejf — ► V%j as K max — > oo and (ii) the eigenvalues, 
Ei(K max ), and eigenvectors of the effective Hamiltonian converge to their limiting values 
faster than the eigenvalues and eigenvectors of the bare Hamiltonian. 

We expressed the HH effective interaction in terms of the matrix element of the iVth 
Jacobi coordinate, i.e. the "last" particle pair, 

<£ Wn,)> = ^^<Wf%A-i)> • (8) 
i<j z 

The relevant hyperspherical degrees of freedom associated with V2 e//(^vi-i) are fjN and 
the hyper angle, 
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r a,a-i , n s 
sin a N = . (9) 

For the construction of V 2 e //(rA,A-i) i n Ref. [0 we made the following ansatz for the 
hyperspherical "2-body" Hamiltonian 

H 2 (p) = ^-^f + V(V2psma N • m) . (10) 
2m p z 

It is a natural choice since Kff contains the canonical kinetic energy associated with the 
two-body variables and f) N (see Eq. (||)). Such an H 2 is in fact an A-body Hamiltonian. 
It contains the hyperspherical part of the A-body kinetic energy and thus becomes a function 
of the collective coordinate p. The hyperradial kinetic energy operator has not been included 
in H 2 . The reason is that we can use a complete basis set for the p-space and therefore we 
do not need to define an effective interaction for the hyperradial part. 

Due to the collective coordinate, p, in H 2 one has automatically a confinement of the 
2-body-system: for moderate values of p the relation < ta,a-i < V%P ensures localization 
of the 2-body wave function and for large values of p the Hamiltonian coincides with the 
bare one, since the NN interaction vanishes. Therefore large overlaps between the model 
space states and the eigenvectors of the 2-body problem are ensured and thus, different from 
the HO effective interaction, there is no necessity for an additional confining potential. 

In order to diagonalize H 2 we have to recouple the last 2 particles into good spin s, 
isospin t, t z and j states 

(T 2 T 3 . . . TaTaz({[Kn-i\Ln-i] £n)K n L n ; S 2 S 3 . . . Sa^JJz\ 

H 2 IT2T3 . . . T'/Ta^K'x^L'x^ £' n )K' n L' n ; S' 2 S' 3 ...S'^jJ J z ) 
= 8[K lf - 1 )[K' fl _ 1 \8s2,s' 2 8s a ,s' 3 ■ ■ ■ 8 Sa . 2 ,s' a _.5t 2 ,t! 2 ■ ■ ■ 5t a „ 2 ,t* a _ 2 



X E E (-1) 2Sa - 2+Sa+S '^(2S a ^ + 1)(2S' a _ 1 + 1)(2s + 1)(2s' + 1) 

ss' j Ja-2 



X 



2 £>A S J 1 2 ^4 S ' 



'(2L N + 1)(2S A + 1)(2J A _ 2 + l)(2j + r 



xJ(2L' N + 1)(2S' A + 1)(2J A _ 2 + l)(2j + 1) 



Ln- 


-1 


£n 


L N 


S A - 


-2 


s 


S A 


Ja- 


-2 


J 


J 


Ln- 


-1 


f 

N 


ji 

N 


Sa- 


-2 


s' 


S' A 


Ja- 


-2 


J 


J 



X 



^ ( _ 1)2 t a _ 2+ t a+ tW (2Ta _ i + i)(2TA_ 1 + l)(2t + 1){2V + 1) 
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-A 



X /, (Ta-2T(A-2)zTaTaz 1 tt z ) (Ta- 2 T(A-2) z T'a Taz 1 t't z ) 

x(K N (£ N ;s)jtt z \ H 2 \K' N (e' N ; s')j t't,)^ (11) 

where S n (T n ) is the total spin (isospin) of the particles 1, 2,..., n. The recoupling leads to 
the following result for the matrix elements of H 2 depending on the A-body HH functions 
and on the spin-isospin states of the last two particles 
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where 



Here 



and 



(K N (£ N ; s)j tt z \ H 2 \K' N (£' N ;s')j = 
x x x x 1 K N (K N + 3N-2) , r K N -i.j , x / 1Q x 



V 



^iSw^(p) = / dn^jyW^sina^)^] , (14) 



= ((In\ s)j tt z \0^\(£' N ; s')j t%) , (15) 



is the usual 2-body matrix element for the operator One sees that H 2 is diagonal in 

the quantum numbers [K N _i], j,t z and 

normally, as in case of the potential models considered in this work, also diagonal in s 
and t. 

Due to the hyperangular integration H 2 explicitly depends on quantum numbers of the 
residual system, i.e. i^jv-i, while it is independent of the other quantum numbers in [Kjv-i]- 
As a result the below defined HH effective interaction depends on the state of the residual 
A — 2 particle subsystem. Such a "medium correction" of the 2-body force is of course a 
great advantage. On the other hand one has to pay for it with greater numerical effort, 
since the effective interaction has to be calculated for all the various states and in addition 
it also depends on the specific A-body system considered. It is similar to the HO multi- 
valued effective interaction [|1^,|13|]. We solve the hyperradial equation on a grid, where H 2 
is diagonalized for each grid point pi and for all the possible values of K N _i,j, s, t, t z in our 
model space. 

Having defined H 2 we employ in Ref. || the Lee- Suzuki [14] similarity transformation 



method to construct the effective interaction. We follow an analogous procedure as that 
of Barrett and Navratil [|I|] for the HO effective interaction. To this end the eigenvectors, 
and eigenvalues, {ej}, of 
H 2 (p), given by Eq. [12|, are used to construct the effective interaction. Denoting by 
\a) the HH functions that belong to our model space, i.e. the HH function |[-?Gv]) with 
Kn < K max , and by \(3) the states that belong to the Q space, Q = {\[K N ]) ; K N > K max } 
the Lee-Suzuki effective Hamiltonian takes the form, 

PH 2 P = PH 2 P + PH 2 QujP , (16) 

where the transformation operator u = QujP is given by the equation 

(/?|i) = 5>Ma)(a|i>. (17) 

a 

If rip is the number of model-space HH basis functions that belong to the subspace [Kjv_i], 
we may solve Eq. flPT|) for uj by choosing a set, A, of rip eigenvectors with the lowest 
eigenvalues |i) and inverting the matrix (a\i). The resulting effective 2-body Hamiltonian 
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n p 



{a\H 2 {K N _ 1 ,p)\a') = Y. 



(a\i)ei(i\a') + ^2(a\i)ei(i\(3)((3\uj\a) 

P 



(18) 



will have the property that P\i), \i) € A, is a right eigenvector of H 2 with eigenvalue e*. The 
effective Hamiltonian is in general a non-hermitian operator, however it can be hermitized, 



using the transformation [|Tj| 

H 2 eff = [P(l + ^uj)P} 1 ' 2 H 2 \P(1 + J^P}- 1 ' 2 . 



(19) 



The effective interaction can be now deduced from H 2e ff, by subtracting the kinetic energy 
term, 



V eff - H 2 eff 



1 K 2 

1 N 

2m p 2 



(20) 



III. NUMERICAL RESULTS 

We apply the EIHH method to the three- and four-body nuclei calculating ground state 
energies (Eb), matter radii (< r 2 > 1 ' 2 ), and two-body correlation functions. We have chosen 
to perform the calculation with two standard noncentral potential models, AV6 and AV14 
0. They are frequently used in other few-body approaches and thus allow to compare to 
other results in the literature. The model AV14 is a fully realistic potential, while AV6 
contains besides spin and isospin dependent central force only the tensor forces in addition. 

In Fig. 1 we show the triton binding energy as function of the hyperspherical quantum 
number K. The convergence pattern looks very similar for both noncentral potential models. 
One sees that quite precise results are already obtained with rather small K values. In fact 
deviations from the converged results amount to less than 0.1 MeV and 0.01 MeV for K > 8 
and K > 14, respectively. It is evident from the figure that the additional potential terms of 
AV14 do not have any important impact on the convergence behavior. In Fig. 1 we illustrate 



in addition our EIHH results from Ref. |2| for the central potential model MTI-III [|16j. The 
comparison with the AV6 result shows that the tensor force does not lead to a deterioration 
of the convergence. 

In Fig. 2 the results for the triton radius are shown. One sees that also in this case 
an extremely good convergence behavior is obtained. For K > 6 one has only very small 
deviations of less than 0.005 fm from the converged AV6 and AV14 results. As for the binding 
energy additional potential terms of AV14 do not affect the convergence behavior. Again 
we show the MTI-III results. Contrary to the case of the binding energy the convergence is 
slowed down a bit by the tensor force. 

Binding energy and radius are typical examples of observables which are not very much 
affected by high momentum components of the wave function. One has to realize that our 
effective interaction approach can be considered in some way as a momentum expansion, 
since with increasing K max higher and higher momentum components can be incorporated in 
the wave function. To study this question better we show in Fig. 3 the two-body correlation 
function 
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p{r) = ($\5{r - r A ,A-i)\^) , 



(21) 



where \I/ is the ground state wave function (/ drp(r) = 1). One sees in Fig. 3a that p(r) 
exhibits a rather broad and asymmetric maximum at about 2 fm. The convergence pattern 
of Fig. 3b shows that at larger distances convergence is obtained with rather low values of 
K. Further increase of K improves the convergence at shorter and shorter distances. In 
Fig. 3b we also illustrate a result with the bare interaction. For an intermediate distance 
between about 1 and 3 fm one finds a slight overestimation of p(r), while at smaller and 
larger r the correlation function is considerably underestimated. 

For the three-nucleon system it is possible to take into account rather high K values 
in the calculation, in fact, e.g., for the highest value we considered (K = 32) one has 
374 hyperangular basis functions. Adding just one further nucleon leads to a considerable 
enhancement of the number of basis functions. Therefore our 4 He calculation is carried out 
only up to K = 18, where one already has 1975 hyperangular basis functions. 

The results for the 4 He binding energy and radius are depicted in Figs. 4 and 5, respec- 
tively. One notes that the convergence behavior is again very good. On the other hand the 
results do not converge as rapidly as for the three-nucleon case. In fact with K = 8 one 
still lacks about 1 MeV binding, while only about one tenth of this value is missing for the 
triton case at the same K. For the convergence of the radius the differences are smaller 
between the 3 H and 4 He results. For K = 6 the 4 He radius is overestimated by 0.015 fm, 
which is only about three times larger than the corresponding overestimation for the triton 
case. However, also for 4 He one finds a rather good convergence for the higher K values. 

From Figs. 4 and 5 it is evident that there are no important differences between the 
convergence patterns of the AV6 and AV14 cases showing again that the additional AV14 
potential terms have a rather unimportant effect on the convergence behavior. Also for 4 He 
we illustrate the results of the central potential model MTI-III in the figures. They exhibit 
a little better convergence behavior than the potential models with tensor force. 

In Fig. 6 we show the two-body correlation function of 4 He. The convergence pattern is 
not as good as for the triton case of Fig. 3. Even at larger distances one needs a higher value 
of K in order to reach the same precision. The lower convergence rate can be attributed to 
the higher density of the alpha particle. 

In Table I we compare our results for the three-nucleon system with other results in 
the literature. Inspecting the table one sees a very good agreement among the various 
methods. It is evident that the binding energy can be calculated with a precision of about 
0.01 MeV in most methods. For the alpha particle we obtain the following results: Eb = 
25.61(4) MeV, < r 2 > 1 ' 2 = 1.4739(4) fm (AV6, no Coulomb force) and E B = 24.34(5) MeV, 
< r 2 > l / 2 = 1.5297(3) fm (AV14, with Coulomb force). The comparison with other results 
in the literature is not as good as for the triton case, e.g., comparing with the SVM result 
2"0fl one finds a difference of 0.2 MeV for the binding energy with AV6. One has similarly 



large differences comparing our Eb for AV14 with that of Ref . |23| . It shows that further 
investigations of the various groups p]j2], ^ , p0| , ^l] , p^ , ^| have to be made in order to reach 
a better agreement. Such a bench mark study is currently under way. 



S 



IV. CONCLUSION 



In this work we have applied the recently developed hyper spherical effective interaction 
method to noncentral potential models including the realistic NN force AV14. In this ap- 
proach the two-body effective interaction depends on the A-body hyperradius and on the 
state of the A-2 rest system explicitly. Our results for the ground state properties of three- 
and four-nucleon systems show that the method leads to an excellent convergence of the 
hyperspherical expansion. The convergence is extremely fast in case of triton, while it is 
a bit slower for the alpha particle. The difference is explained by the higher density of 
the four-nucleon system, since our approach can be interpreted as a kind of momentum 
expansion. 

In view of further applications it will be important to incorporate also three-nucleon 
forces in the EIHH method. Therefore we would like to mention that, similar to the HO 
approach |]T|,PE| , the present formalism can be extended to derive an HH three-body effective 
interaction. 
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FIGURES 



FIG. 1. Triton binding energy as function of the hyperangular quantum number K with various 
NN potentials models: AV14 (a), AV6 (b), and MTI-III (c). The value for the largest K is indicated 
by a dashed line. 

FIG. 2. Triton matter radius as function of the hyperangular quantum number K with various 
NN potentials models with notations as in Fig. 1. 

FIG. 3. EIHH results for the two-body correlation function of triton with the AV14 potential 
with various K relative to the K = 26 result; also shown is the K = 26 result with the bare 
interaction. 

FIG. 4. Alpha particle binding energy as function of the hyperangular quantum number K 
with various NN potentials model with notations as in Fig. 1. 

FIG. 5. Alpha particle matter radius as function of the hyperangular quantum number K with 
various NN potentials model with notations as in Fig. 1. 

FIG. 6. EIHH results for the two-body correlation function of 4 He with the AV14 potential 
with various K relative to the K = 16 result. 
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TABLES 



TABLE I. Comparison of binding energies (Eb) in [MeV] and root mean square radii 
(< r 2 > x / 2 ) in [fm] obtained with the present effective interaction method in the HH formal- 
ism (EIHH) with results of other methods. For EIHH the number in parenthesis indicates the 
variance with respect to the result obtained with K = K max — 2. The quality of the convergence 
can be inferred from Figs. 1, 2 







AV6 




AV14 




Nucleus 


Method [Ref.] 


E B 


< r 2 >l/2 


E B 


< r 2 >l/2 


3 H 


EIHH 

Faddeev Q 
GFMC || 
VMC |19| 
SVM |20| 
HH [21] 
CHH H 


7.1602(3) 

7.15 
7.22(0.12) 

7.15 


1.77475(5) 
1.75(0.10) 
1.76 


7.6814(3) 
7.670 

7.53 

7.684 
7.69 


1.77615(3) 
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1.79 



AV14 




